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ABSTRACT 

The rotationally resolved spectrum of the <— A X E+ 000-000 transition of C3, centered at 

4051.6 A, has been observed along 10 translucent lines of sight. To interpret these spectra, a new 
method for the determination of column densities and analysis of excitation profiles involving the sim- 
ulation and fitting of observed spectra has been developed. The populations of lower rotational levels 
(J<14) in C3 are best fit by thermal distributions that are consistent with the kinetic temperatures 
determined from the excitation profile of C2. Just as in the case of C2, higher rotational levels (J>14) 
of C3 show increased nonthermal population distributions in clouds which have been determined to 
have total gas densities below ^500 cm~ 3 . 

Subject headings: astrochemistry — ISM: lines and bands — ISM:molecules 
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1. INTRODUCTION 

The existence of carbon chain molecules in diffuse in- 
terstellar clouds raises questions about the chemistry 
that occurs in these clouds as well as the fate of large car- 
bonaceous molecules produced throughout the life-cycles 
of stars. Carbon chain molecules have been suggested as 
the source of t he confounding diffuse interstellar bands 
({Douglas) 119771) . although efforts to identify a specific 
molecular carrier have so far failed. C7 and I-C3H7 have 
been the most recently rejected (|Mc£jalle^dJ [2002b, 
120011) in a long list of candidates |l Icrhiij jlj.^) . The 
smallest carbon chain molecule, C2, has been used as a 
powerf ul probe of the temperature a nd density in diffuse 
clouds l)van Dishoeck fc Blacldll"989l) and, until very re- 
cently, diffuse cloud chemistry was understood only in 
terms of atoms and diatomics. 

After repeated efforts llClegg fc Lambertl Il982t 
lHaffner fc Meverlll995l: iSnow et alJll988|) . the first clear 
det ections of C 3 in th ree diffuse clouds were reported 
by iMaier et alJ pOOlfr . Since this first observational 
study of C3 in the diffuse interstellar medium, a 
detailed m odel has been deve loped to understand its 
excitation ijRoueff et al.l 12002)) and a low resolution 
survey has dete cted it in 15 translucent sight lines 
llOka e t al. 2003). Along with infrared measurements 

nmiir 



die et al] 119991 iMcCall et alJl!998l l2002alh 

which can dire ctly probe the inter stellar cosmic ray 
ionization rate i jMcCall et al.ll2003bD . the detections of 
C3 demonstrate that the chemistry of diffuse clouds may 
be significantly richer than expected. 

The "A4050 gro up" of C^ was first observed in emission 
in comet Tebutt l)HuggmsHl881)) . and has sinc e been ob- 
serve d in the photospheres of many co ol stars lUorgensen 
1994), in a few protoplanetary nebulae ( Hrivnak & Kwok 
1999), and toward many other comets (see, for exam- 
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pie, recent work bv lRousselot et alJ 1)20011) 1. The A4050 
group was fir st observed in the laboratory by C. W. 
Raffctv (1916) i n emis sion from the flame of a Meeker 
burner. iDouglasI l) 195 It) identified the carrier of the group 
as C3, and a detailed assignment as A 1 ^ «— A 1 !]^ was 
made bv lGausset et al.l l)1965|) . At longer wavelengths, 
C3 has also been observed in the material around the 
asymptotic giant branch star IRC +10216 using both 
the asymmetric vi brational stretching mode at A.9fim 
(Hink le et al.lll988|) and the anomalously low fre quency 
bending mode at 157/um ijCernicharo et al.ll2000)h 

As a linear carbon chain, like C2, C3 has no perma- 
nent dipole moment and thus lacks the ability to cool 
efficiently via rotational transitions. High rotational lev- 
els of C3, once produced during chemical formation or 
following UV pumping, can only relax slowly by colli- 
sional de-excitation an d have significant l ifetimes under 
interstellar conditions. IMaier et al.l l)2001[) reported that 
the low J levels (J<14) in their spectra of C3 toward £ 
Oph were fit by a Boltzmann distribution at 7) =60 K, 
whereas J>14 were fit by a Boltzmann distribution at 
7^=230 K. 

Depending on the intensity of the UV radiation field, 
the rates of collisional (de-)excitation an d UV pump- 
ing, a nd the photochemical lifetime of C3, Roucff et al. 
( 2002) suggest that the high rotational levels of C3 may 
retain information about the temperature at which C3 
was formed as well as the tem perature and density of 
the cloud. iRoueff et al.l l)2002fl applied their sophisti- 
cat ed C3 excita t ion m odel to all three of the detections 
by IMaier et alJ (|2001|) - £ Oph, £ Per, and 20 Aql - 
along with their own observations of HD 210121, but 
none of the observations determined the population of 
C3 in J > 14 well enough to constr ain the full excitation 
model. iGalazutdinov et "all 1)2002)) have since provided 
spectra of C3 toward \ Oph and HD 152236, without 
analysis of the excitation. 

With the guidanc e of the low resolution survey results 
of lOka et alJ l|2 p03[[) we have measured the spectra of 24 
reddened stars at high resolution in order to study the 
excitation profile of C3 in a large sample of diffuse clouds. 
Here we present the rotationally resolved spectra of C3 
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Table 1. Data For Program Stars 



Star 


Name 


Type 


V 


E(B-V) 


D (pc) 3 - 


S/N b 


Date (UT) 




e Cas 


B3 III 


3.38 


0.05 


136 


2400 


9009 Dec 1 ^ 


HD 20041 




AO la 


5.79 


0.72 


440 


1000 


2002 Dec 14 


HD 21389 




AO lac 


4.54 


0.57 


940 


1300 


2002 Dec 14 


HD 21483 




B3 III 


7.06 


0.56 


440 


1000 


2002 Dec 15 


HD 22951 


40 Per 


B0.5 V 


4.97 


0.27 


283 


1700 


2002 Dec 15 


HD 23180 


o Per 


Bl III 


3.83 


0.31 


280 


1850 


2002 Dec 14 


HD 24398 


C Per 


Bl lb 


2.85 


0.31 


301 


2200 


2002 Dec 15 


HD 24534 


X Per 


09.5pe 


6.10 


0.59 


590 


1200 


2002 Dec 15 


HD 24912 


£ Per 


07e 


4.04 


0.33 


470 


2600 


2002 Dec 14 


HD 34078 


AE Aur 


09.5 Ve 


5.96 


0.52 


620 


1200 


2002 Dec 14 


HD 35149 


23 Ori 


Bl V 


5.00 


0.11 


450 


2100 


2002 Dec 15 


HD 36371 


X Aur 


B5I ab 


4.76 


0.43 


880 


2000 


2002 Dec 15 


HD 37022 


9 1 Ori C 


06 


5.13 


0.34 


450 


2000 


2002 Dec 15 


HD 37061 




Bl V 


6.83 


0.52 


580 


1300 


2002 Dec 15 


HD 41117 


X 2 Ori 


B2 Iae 


4.63 


0.45 


1000 


2100 


2002 Dec 15 


HD 42087 


3 Gem 


B.25 Ibe 


5.75 


0.36 


1200 


1600 


2002 Dec 14 


HD 43384 


9 Gem 


B3 lb 


6.25 


0.58 


1100 


1600 


2002 Dec 14 


HD 53367 




B0 IVe 


6.96 


0.74 


780 


1300 


2002 Dec 15 


HD 62542 




B5 V 


8.04 


0.35 


246 


700 


2002 Dec 14 


HD 169454 




B1.5 la 


6.61 


1.12 


930 


450 


2002 Jul 14 


HD 204827 




B0 V 


7.94 


1.11 


600 


650 


2002 Dec 14, 2002 Jul 13 


HD 206267 




06 f 


5.62 


0.53 


1000 


1200 


2002 Jul 12, 2002 Jul 13 


HD 207198 




09 Be 


5.95 


0.62 


1000 


700 


2002 Jul 12, 2002 Jul 13 


HD 210839 


A Cep 


06 If 


5.04 


0.57 


505 


700 


2002 Jul 13 



a The distance to 9 1 Ori C is from Shuping and Snow (1997). All 
other distances are taken from Thorburn et al. (2003) and Oka ct 
al. (2003). 

b Signal to noise per pixel, measured near 4049 A. 



at high signal to noise (S/N) in 10 of these translucent 
sight lines. We present a new spectral fitting method for 
determining the populations in each J level and examine 
the relationship between the rotational excitation of C3 
and C2. 

2. OBSERVATIONS 

Observations were carried out using the Shan e 3-m 
telesc ope and the Hamilton Echelle Spectrograph l|VogtJ 
|1987|) at Lick Observatory on 2002 July 12 through 
14 and us ing the High Res olution Echelle Spectrometer 
(HIRES) <)Vogt et al.lll994H on the 10-m Keck I telescope 
atop Mauna Kea on 2002 December 14 and 15. 

The Hamilton spectrograph operates with an echelle 
grating and two cross-dispersing prisms, is mounted at 
coude focus, and provides nearly complete spectral cov- 
erage from 3,500-10,000 A. The detector used was the 
Lick3 L4-9-00AR thinned 2048x2048 CCD with 15^m 
pixels. The l."2 wide and 2."0 long slit provides a re- 
solving power i?=A/AA~60,000, where the resolution AA 
refers to the full width at half maximum (FWHM) of the 
instrument profile, ~0.1 A at 6,000 A, which spans two 
pixels on the detector. The reduced data have a disper- 
sion d=0.0336 A pixel" 1 at 4050 A. 

The HIRES instrument is also a cross-dispersed echelle 
spectrograph, and features a backside-illuminated Tek 
2048x2048 CCD detector functioning over the 3,500- 
10,000 A range and capturing a spectral span of 1200- 
2500 A per exposure. To minimize total readout time for 
our relatively bright targets, high gain mode and dual- 
amp readout were used. For these observations spectra 
in the 3985-6420 A range were recorded with each expo- 
sure, with small gaps in the spectral coverage starting 
at ^5000 A and increasing to ~30 A at 6315 A, due to 



echelle orders falling off the edge of the CCD. A slit size 
of 0."574 x 14" was used to achieve a resolving power of 
i?~67,000. The detector position was constant through- 
out each night and moved only slightly between the two 
nights of observations. The reduced data have a disper- 
sion d=0.0287 A pixel" 1 at 4050 A. 

The data were reduced using standard procedures in 
the echelle package in NOAO's Image Reduction and 
Analysis Facility (IRAF). As many as 120 flat field im- 
ages were taken before and 180 after the observations to 
ensure that the S/N of the reduced spectra was photon- 
counting limited. Calibration and target images were 
bias corrected and combined, before target images were 
flat fielded and echelle apertures extracted. Wavelength 
calibration was performed with standard routines using 
spectra of a ThAr lamp, and one-dimensional spectra 
of the targets were extracted for further analysis. Gen- 
erally, the entire aperture blaze was fit with a low or- 
der polynomial to normalize the spectra. In cases where 
the spectral features of interest were observed very close 
to the center (maximum) of the aperture blaze, only a 
small spectral region around the features was continuum 
fit with a low order polynomial. The program stars and 
their properties, along with the S/N achieved at 4049 A, 
are shown in Table Figure 1 presents the C3 spectra 
observed along the 10 lines of sight with detectable C3 
absorption. 

3. SPECTRUM OF C3 

Unlike molecules which can radiate away excess en- 
ergy, the distribution of population among rotational 
levels for molecules in the ISM without a permanent 
dipole moment, such as C3, is a delicate balance be- 
tween collisional and radiative processes. This has 



ROTATIONALLY RESOLVED C 3 



3 



Table 2. Molecular Data for the A 1 H u <— transition of C3 



Wavelength (A) Line fj>,j» (xl0 3 ) a 



4049.770 


R(24) 


4.24 


4049.784 


R(26) 


4.23 


4049.784 


R(22) 


4.26 


4049.810 


R(20) 


4.30 


4049.810 


R(28) 


4.21 


4049.861 


R(18) 


4.33 


4049.861 


R(30) 


4.20 


4049.963 


R(16) 


4.36 


4050.081 


R(14) 


4.42 


4050.206 


R(12) 


4.48 


4050.337 


R(10) 


4.57 


4050.495 


R(8) 


4.70 


4050.670 


R(6) 


4.92 


4050.866 


R(4) 


5.34 


4051.069 


R(2) 


6.40 


4051. 250 b 


R(0) 


16.00 


(4051.309) 




4051.461 


Q(2) 


8.00 


4051.519 


Q(4) 


8.00 


4051.590 


Q(6) 


8.00 


4051.681 


Q(8) 


8.00 


4051. 741 c 


P(2) 


1.60 


4051.793 


Q(10) 


8.00 


4051.929 


Q(12) 


8.00 


4052.062 


P(4) 


2.66 


4052.089 


Q(14) 


8.00 


4052.271 


Q(16) 


8.00 


4052.424 


P(6) 


3.08 


4052.473 


Q(18) 


8.00 


4052.698 


Q(20) 


8.00 


4052.792 


P(8) 


3.29 


4052.940 


Q(22) 


8.00 


4053.180 


P(10 


3.43 


4053.208 


Q(24) 


8.00 


4053.490 


Q(26) 


8.00 


4053.591 


P(12) 


3.52 


4053.794 


Q(28) 


8.00 


4054.020 


P(14) 


3.58 


4054.112 


Q(30) 


8.00 


4054.458 


P 16) 


3.64 


4054.908 


P(18) 


3.67 


4055.373 


P(20) 


3.70 


4055.877 


P(22) 


3.73 


4056.410 


P(24) 


3.75 


4056.961 


P(26) 


3.77 


4057.531 


P(28) 


3.78 


4058.121 


P(30) 


3.80 



a Here / =0.0160 was used, consistent with |Mai«]^^aIJ_J20pj|) 
andlOka et al.l E00^ ) : IRoueff et alj O002T) and lGalazutdinov et all 
12002ft use f =0.0146. 

b The R(0) assignment by Gaussct et al. ( 1965f) (shown paren- 
thetically) is inconsistent with our observations. We use the ob- 
served value, confirmed by laboratory measurements IMcCall et alj 
2003a}), in the rotational level popul ation model. 

c Calculated transition wavelength (Roueff et al] |2002() . 



been studied in det a il for the homonuclear diatomics H2 
(Blac k~ Dalgarnol II1977D. and references therein) and 
C9 ijvan lDishoe^!rfcBlacM ri982). and more recently a 
simi lar method of ana lysis has been used to examine 
C3 ijRoueff et alJ 12002). In general, these models show 
that the populations of lower rotational levels (J<1A for 
C3) are controlled mostly by the kinetic temperature of 
the gas along the line of sight, whereas the populations 
in higher rotational levels are determined by the com- 
petition between radiative pumping and collisional de- 
excitation. The net result can be a non-thermal enhance- 
ment of high J populations in low density environments 



as observed bv lMaier"et al. (200]]) toward £ Oph. 

In order to determine the total column density of C3 
and measure the population distribution among rota- 
tional levels, we have used two models to fit 40 rovibronic 
transitions between 4049-4055 A. A thermal excitation 
model incorporating either one or two temperatures was 
used, as well as a model that fit the population in each 
rotational level to the observed spectra. These two meth- 
ods provided a means to extract information from over- 
lapping or unresolved transitions in the C3 Q-branch and 
i?-branch bandhead. 

Laboratory measurements l)Causset et al.l Il965j) of 
transition wavelengths were used in both models, with 
the exception of a shifted R(Q) line that was used in the 
rotational level population model (see Section l3~4"|) . 

3.1. Thermal Excitation Model 

In this model the relative population of each rotational 
level ( J<30) is determined using a thermal distribution 
at temperature, T. The fraction of total molecules in a 
particular J level is given by 

where q r is the rotational partition function given by 

00 

q r = ( 2 J + l)e- hcB ^ J+1 ^ kT 

J even 

For C3 the ground state rotational constant Bo=Q. 43057 
cm -1 ( Sc hmuttenmaer et alJllffflfil) . The two parameters 
T and AT(C3), along with the /-values shown in Table 
121 then determine the equivalent width, Wj' t j", of each 
transition using the standard relationship 

Wj>,j» = 8.853 x 10- 21 \ 2 Njf j, j„ A (2) 
where fj',j", the oscillator strength for a given tran- 
sition J' <— J" at wavelength A in A, is given by 
the product of the electronic oscillator strength /„;, 
the Franck-Condon factor, and the Honl-London fac- 
tor. The units for Nj and Wj are cm -2 an d A, respec- 
tively. Here we use /=0.016, consistent with lMaier et alJ 
l)2001|) and lOka et ail l)2003|) . for the product / of the 
electronic o scillator strength a nd Franck-Condon fac- 
tor, w hereas IRoueff et al.l l)2002D and lGalazutdinov et alJ 
(2002) use /=0.0146. Future determinations of the os- 
cillator strength and Franck-Condon factor will scale the 
derived column densities accordingly. The Honl-London 
factors for this band (A = 1 <— 0) arc (J + 2)/2(2J + 1), 
1/2, and (J-1)/2(2J+ 1) for the R-, Q-, and P-branch 
transitions, respectively. Table shows the adopted 
wavelengths, assignments and fj\j» for each transition. 

Beginning with an approximate estimate of the 
C3 excitation temperature (T~45 K), column density 
(A^(C 3 )-10 12 cm" 2 ), and line width (FWHM=5.2 km 
s _1 ) a simulated spectrum was produced. Standard 
curve-fitting routines were then used to minimize the dif- 
ference between the observed and simulated spectra, by 
varying these three parameters (T, N(Cs), and FWHM). 
If a second temperature component was observed, in the 
form of a pronounced i?-branch bandhead that was not 
reproduced with the single temperature simulation, then 
a spectrum incorporating two temperatures XJ and T^, 
and column densities Ni(Cs) and Nh(Cs) was used. 
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Fig. 1. — Spectra of translucent sight lines toward which C3 was detected, along with thermal excitation model fits. Note how the 
Q-branch at 4051.6 A is very well fit, while the R bandhead at 4049.8 A is sometimes underestimated. Spectra are all shifted to rest 
wavelengths and the K I line at 4047.2 A is marked by a vertical line. Broad stellar lines have been removed from the spectra of ( Per, HD 
21483, AE Aur, HD 36371, and HD 53367 using Gaussian line profiles. 



Table 3. Best Fit Parameters for Thermal Excitation Model of the 
Observed Spectra 



Star 


JV(C 3 ) a 


^((C 3 ) a 




Ti 


T h 


FWHM b 
(km s — 1 ) 


HD 21483 


4.60 ± 0.28 


2.60 


2.00 


28.3 ± 1.3 


150 ± 13 


5.6 ± 0.1 


HD 23180 


1.27 ± 0.13 






127.3 ± 3.9 




5.2 ± 0.2 


HD 24398 


1.30 ± 0.35 






132 ± 10.3 




5.1 ± 0.4 


HD 24534 


1.77 ± 0.19 


0.49 


1.28 


40 ± 6.1 


250 ± 41 


3.7 ± 1.0 


HD 34078 


6.21 ± 0.68 






301.4 ± 15.4 




7.2 ± 0.3 


HD 36371 


0.75 ± 0.06 


0.13 


0.62 


20.1 ± 2 


250 ± 52 


3.1 ± 1.2 


HD 53367 


2.08 ± 0.26 






60.8 ± 1.5 




5.6 ± 0.1 


HD 62542 


10.37 ± 0.53 


7.35 


3.01 


75.5 ± 7.2 


200 ± 76 


5.6 ± 0.1 


HD 169454 


2.24 ± 0.66 






23.4 ± 1.4 




6.5 ± 0.4 


HD 204827 


11.13 ± 0.87 






40.6 ± 0.8 




8.5 ± 0.2 



a Column density N in 10 cm . 

b Full-width at half maximum derived from direct fit to the ob- 
served spectrum; instrumental line width at 4050 A is ~4.5 km 
s _1 . 



The initial estimate for the line width was larger than 
the instrument resolution (~4.5 km s _1 ) and character- 
istic of the line widths in the observed spectra (Table 
|3J). Broad lines are observed because translucent sight 



lines typically sample multiple diffuse clouds, which can 
have velocity separations co mparable to, or larger than, 
our instrumental line width l|Weltv fc Hobbsfe OOl). The 
velocity distribution of the clouds, convolved with the 
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instrumental line width, produces observed features that 
in some cases (e.g. HD 204827) are rather wide. Since we 
are observing through multiple clouds, the temperatures 
derived from our model represent an integrated average 
over all the clouds along the line of sight. In the opposite 
extreme, the best fit line widths for X Per and x Aur are 
slightly below the instrument profile due to the relatively 
large level of noise in those spectra. 

Figure ^ shows the observed spectra and best fit ther- 
mal excitation model for 10 detections, with model re- 
sults summarized in Tabled Measurement uncertainties 
are discussed in Section l3~3l 

3.2. Rotational Level Population Model 

The most common method for determining the pop- 
ulation of a particular rotational level J is to measure 
the equivalent width W of an unblended absorption fea- 
ture. If multiple transitions (A J = 0, ±1) from the same 
lower J state are observed, then each transition probes 
the same population and multiple measurements of the 
same quantity may be made. 

In the case of our C3 spectra, the low J transi- 
tions (J<12) in the i?-branch are well resolved, but the 
stronger Q-branch lines starting from the same J lev- 
els pile atop one another. Similarly, transitions from 
higher J levels that are well resolved in the Q-branch 
(e.g. Q(16) and Q(22)) are blended into a single fea- 
ture in the i?-branch bandhead. In order to extract in- 
formation from all of the observed transitions, a model 
was used which varies the population in each J level as 
a free parameter, produces a simulated spectrum, and 
then minimizes the difference between the observed and 
simulated spectrum. 

The initial population values for the model were de- 
termined from the best fit thermal distribution. A curve 
fitting routine was then used to vary each Nj and min- 
imize the residuals. After fitting the population in each 
J level, a Boltzmann plot was generated for each source, 
and kinetic temperatures which best fit J<14 and J>14 
were determined. Figure |3 shows the observed and mod- 
eled spectra, with Boltzmann plots derived from the fit. 
The detection limits for each J are shown as dotted lines 
in the Boltzmann plots. Results from the model are listed 
in Table 21 with uncertainties discussed below. 

3.3. Uncertainties and Detection Limits 

Uncertainties in the measurement of column densities 
for individual rotational levels are dependent on the un- 
certainty in the measurement of the equivalent width W . 
The la uncertainty in W is given by 

SWio = (wdf^a (3) 

where w is the FWHM of the observed transition (in A) , 
d is the dispersion of the instrument (in A pixel -1 ), and 
cr is the standard deviation of the normalized continuum- 
that is, (S/N) -1 . The 5W for each transition is converted 
to an uncertainty in population, SNj using Equation 2. 
For a given J the transition with the largest oscillator 
strength will give the smallest uncertainty in Nj. In all 
cases except the R(0) transition, the Q-branch has the 
largest oscillator strength for a given J (see Tabled and 
this is the /-value we use for our uncertainty estimate. 



The detection limits for the total column of C3 are 
dependent on the excitation of C3. A larger total popu- 
lation of C3 may go undetected if it is distributed among 
many rotational levels than if the population is all con- 
fined to a few low J rotational levels. For our upper 
limits listed in Table El we assume that molecules are all 
distributed evenly among J=0 - 20. Although the choice 
of 1 1 levels among which the distribution is evenly spread 
is not based on any underlying physics, a comparison 
with the spectrum of £ Per shows that for the purposes 
of estimating our detection limit, such an approach is not 
unreasonable. 

3.4. The R(0) Transition 

During the initial analysis of the Boltzmann plots it 
became evident that R(0) was uniformly underpopulated 
relative to what was expected from a thermal distribu- 
tion that fit J=2 - 10. A close inspection of the spectra 
revealed that the observed R(0) transition was signifi- 
cantly shift ed toward the blue f rom the laboratory as- 
signment of iGausset et all £1965). This was particularly 
apparent in the spectrum of HD 62542, and can be noted 
by comparing the t herm al excitation model (which uses 
the IGausset et all l)196jJl assignment of R(Q)) with our 
observations in Figure [D Analysis of th e higher resolu- 
tion spectra of C Op h (Maie r et al]l2001fl and HD 152236 
l)Galazutdinov et alJl2002f) shows similar indications of a 
spectral shift, however, neither dataset is of sufficient 
S /N to conclude that there is a significant difference be- 
tween the observed spectrum and the laboratory assign- 
ment. Motivated by this discrepancy, a group including 
two of us (MA and BJM) revisited the laboratory spec- 
trum of the 4051 A b and of C3 with high er resolution 
than our observations ijMcCall et alJl2003aj) . Indeed, we 
found the R( 0) transition t o be b lue-shifted from the 
assignment of IGausset et al.l l)1965T and consistent with 
our astronomical spectra. Since R(0) and P(2) share the 
same upper state, our measurement of the R(0) transi- 
tion implies that the calculated position of the P(2) line 
(Table 0) will shift to the blue and overla p Q(8). This 
is confi rmed by laboratory measurements ijMcCall et alJ 
I2003aft . 

4. RESULTS AND ANALYSIS 

4.1. Column Density Comparisons 

Table |5l summarizes our measurements of N(Cs) along 
with previous work. In their determination of C3 col- 
umn densi t ies fr om an extensive, low resolution survey, 
lOka et al] l)2003j) measured the equivalent width of the 
unresolved Q-branch and assumed that it contained half 
the total intensity of the C3 band. While this is a good 
approximation for a thermal distribution at T<50 K, 
we note that this assumption underestimates the total 
column density when there is significant population in 
higher J levels. Comparisons of N(C 3 ) for HD 204827 
and HD 169454 show either agreement w ithin error, or a 
slight overestimate bv lOka et al J l)2003[) . Notably, these 
are the two clouds for which we measure only a low tem- 
perature component in the excitation profile. On the 
other hand, AE Aur and HD 21483, which show pro- 
nounced . R-branch bandhea ds (Figure GJ, are underesti- 
mated bv lOka et alJ l)2003|) due to a significant popula- 
tion of C3 in high rotational levels. 
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Fig. 2. — Spectra of C3 detections with rotational excitation model fit (offset for clarity) and the Boltzmann plot (left) derived from fit 
to individual rotational levels. The dotted curve depicts the lc detection limit as a function of J. Points without error bars indicate upper 
limits. Dashed lines indicate fits to datapoints for J<14, and J>14, where the slope is inversely proportional to kinetic temperature. Note 
that the spectra are on different intensity scales and the very large columns of C3 toward HD 62542 and HD 204827. 



Our determinat ion of N(G^) for C P er is larger than the 
measurements of IMaier et alJ ((2001) due to our higher 
S/N, which allows us to measure the populations in 
higher J le vels. The col u mn d ensity determination for 
C Qph bv IRoueff et alJ IpOO^ . using the data from 
Maicr et al.l <|2001|> . is consistent with our measurements 
and highlights the utility of their model for determining 
N(C 3 ). 



4.2. Correlations 

The correla t ion b etween C2 and C3 pointed out 
by lOka et all l|200^ may be tested more rigorously 
with the precise determinations of the C3 column den- 
sity presented here. The N(C2) measurements of 
iThorburn etaD (2003) are plotted against our measure- 
ments of AT(C3) in Figure El 

The most notable outlier is HD 62542, which has an 
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Fig. 2. — (continued) 



extremely high column density of C3, yet an average to 
low column of C2. The intriguing line of sight toward 
this star seems to be a highly unorthodox cloud; not 
only does it seem to b e "missing" C2, it has an extreme 
UV extinction curve (Cardclli fc Savage! 119881) . a n un- 
usually low column of CH + l)Cardelli et al.l 11990). and 
it lacks the diffu se bands (DIBs) characteristic of red- 
dened sightlines llSnow et alJl2002f). Following the anal- 
ysis of lCardelli et all ljl99(1) . ISnow et all l|2002ft have in- 
terpreted the lack of DIBs toward this star by suggesting 
the sight line is dominated by a dense cloud stripped of 
its diffuse outer layers. Considering the weak DIBs and 



the seemingly low column of C2, it is noteworthy that 
a correlation between C2 and some DIB s, such as A4963 
and A 5769, has recently been identified ijThorburn et alJ 
12001 . However, the "C 2 DIBs" are in general much 
weaker than the classical DIBs, such as A6284 and A5797, 
which have been reported to be unusually weak toward 
HD 62542 <|Snow et alJl2002|) . Nonetheless, one might 
speculate that the processes which stripped the cloud of 
its strong DIBs may also have destroyed some of its C2. 
On the other hand, it is unclear how C3 would survive 
a process that destroys C2, so one could suggest instead 
that the conditions in this cloud may favor carbon chain 
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Table 4. Rotational Level Populations Determined by Fitting the 
Observed Spectra 



NjixlO 11 cm" 2 ) 



./ 


HD 21483 


HD 23180 


HD 24398 


HD 24534 


HD 34078 


HD 36371 


HD 53367 


HD 62542 


HD 169454 


HD 2048: 





1.22 


0.09 


0.34 


0.14 


0.88 


0.15 


0.94 


2.64 


0.49 


1.83 


2 


5.79 


0.92 


1.01 


1.04 


2.55 


0.27 


2.83 


6.46 


2.13 


15.42 


i 


7.67 


1.21 


1.54 


1.35 


3.65 


0.42 


3.04 


10.04 


4.66 


22.54 


D 


O.o4 


1 1 Q 

1.18 


i fin 


O 1 A 

z. 14 


A f^l 
4.01 


n a o 
U.4z 


Q OQ 

o.yo 


iz.yo 


4. i o 


24.2 ( 


8 


6.19 


1.26 


1.84 


1.72 


7.06 


0.21 


3.17 


15.15 


2.27 


19.44 


10 


3.87 


1.43 


1.61 


1.49 


7.53 


0.29 


2.34 


11.59 


< 0.96 


11.68 


12 


2.61 


1.29 


1.84 


1.28 


8.03 


0.18 


2.17 


11.04 


< 0.96 


6.53 


14 


1.43 


0.82 


1.30 


1.36 


7.85 


< 0.15 


0.57 


9.22 


< 0.96 


4.24 


16 


1.42 


1.19 


1.44 


1.31 


5.30 


0.17 


0.67 


6.94 


< 0.96 


3.48 


18 


1.90 


0.82 


1.39 


1.26 


3.57 


0.49 


< 0.31 


5.54 


< 0.96 


< 0.76 


20 


1.26 


0.82 


0.91 


1.10 


2.38 


0.32 


0.51 


4.36 


< 0.96 


< 0.76 


22 


0.90 


0.73 


0.78 


0.87 


3.65 


0.26 


< 0.31 


3.73 


< 0.96 


< 0.76 


2d 


1.43 


0.82 


0.99 


0.97 


2.74 


0.19 


< 0.31 


2.74 


< 0.96 


< 0.76 


26 


0.91 


1.13 


0.62 


0.76 


1.83 


< 0.15 


< 0.31 


1.19 


< 0.96 


< 0.76 


28 


1.14 


0.72 


0.44 


0.91 


2.28 


< 0.15 


< 0.31 


0.73 


< 0.96 


< 0.76 


30 


1.28 


0.64 


0.52 


1.12 


1.83 


0.44 


0.42 


0.64 


< 0.96 


< 1.15 


SNj* 


0.40 


0.21 


0.17 


0.27 


0.38 


0.15 


0.31 


0.57 


0.96 


0.76 



a The l<r uncertainty SNj is determined using Eqns 2 and 3 and 
/=0.016. Note that SNj for R(0) is half of the listed value. 



formation. Due to the peculiarity of the line of sight to- 
ward this star, we have excluded it from the correlation 
analysis. 

Fitting all of our targets except HD 62542 shows 
a correlat ion between C2 and C3 similar to that ob- 
served by O ka et alJ (|2p03). The correlation coeffic ient 
for this work is 0.898, whereas lOka eta!] lj200l re- 
port 0.932. HD 169454 appears about 2a above the 
tre ndline in Figure IS) but the value of N(Cq) reported 
by iThorburn et alJ (120031) is significantly larger than 
the values reported byly an Dishocck fc Black! l|1989t) and 
IGredel MvmThl ill 9861) . which are more consistent with 
the correlation measured here. AE Aur falls significantly 
below the trendline; perhaps the same processes that 
could have destroyed the DIBs and C2 in HD 62542 
have occurred, to a lesser extent, in AE Aur l(Snow et al.l 
2002). Another consideration is the variability with 
time of the molecular species toward AE Aur, where 
iV(CH) has increase d by 20% over the last 10 years 
(Rolli nde et al.ll2003|) . However, there were only 2 years 
between the observations of C2 and C3 reported here, 
while either N(C2) is a factor of ^2 below or N(Cs) 
is a factor of ^2 too large relative to the general trend 
observed for the other sight lines. The possible enhance- 
ment of N(Cs) is significantly larger than expected from 
the time variability of A(CH) alone. 

Since HD 204827 has such a large column of C3 rel- 
ative to the rest of the sample, much of the correlation 
seems to rest on this datapoint. With the high sensitiv- 
ity measurements of iV(Cs) presented in this work, the 
correlation between C2 and C3 toward stars with rela- 
tively low column densities (Af(C3)<4xl0 12 cm -2 ) may 
be tested. Excluding sources with A r (C3)>4xl0 1 W 2 
(i.e. HD 204827, HD 62542, AE Aur, and HD 21483) 
and HD 169454 for reasons discussed above, we find a 
correlation coefficient of 0.841, supporting the correla- 
tion determined for the entire sample. 

The observed correlation between C2 and C3 is strik- 
ingly similar to the correlation between C2H and C3H2 



measured by iLucas fc Lisztl l)2000D . In many cases the 
ratio of N{C 2 )/N(C 3 ) listed in Table falls within the 
mean ab undance ratio ( C2H /0-C3H2 } = 27.7±8 mea- 
sured bv ILucas & Lisztl l|2000D . It is unclear if this is co- 
incidence or if there is some underlying chemistry which 
links these species and their relative abundances. 

4.3. The Excitation of C2 

The excitation profile of C3 was compared to the phys- 
ical conditions determined by the well-understood rota- 
tional excitation of C2 . Table El summarizes previous 
observations in the literature of rotationally resolved C2 
toward the stars where we detected G3. To the exist- 
ing measurements we have added our C2 spectrum of 
HD 204827 measured at Lick and the spectrum of X Per 
measured at Apache Point Observatory (Thorburn 2003, 
private communication) . 

All measured equivalent widths were converted to pop- 
ulations using the standard relationship in Equation 2. 
For multiple measurements of the same J level, weighted 
averages were used for the population. The rotational 
populations were then simulated using t h e full C2 ex- 
citation model of Ivan Dishoeck & Black! l|1982|) . prob- 
ing the density-temperature parameter space that deter- 
mines the excitation of C2 while minimizing the differ- 
ence between the model and observed values 4 . In this 
way a kinetic temperature T(C2) and a derived density 
of collision partn ers, n r .oii , was determined for a ll targets. 

As detailed by van Dishocck & Black (1982), the de- 
rived density n co u is dependent on the collision cross- 
section Co °f C2 with H2 or H and a scaling factor I 
for the standard radiation field, which we assume to be 
2xl0 -16 cm 2 and 1, respectively. Changes in these val- 
ues will scale n co ii accordingly. 

Figure 0] ( left) shows a comparison of the temperature 
Tiow(C3) determined by fitting the low (J<14) rotational 

4 A web-based C2 calculator provided by BJM is available at 
|http: //dibdata. org/| 
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Table 5. Observed Properties of C3 and C2 



Star 


Name 


N(C 3 ) 




Thigh (C3) b 


N(C 2 ) C 


T(C 2 ) d 


n C oii(C 2 ) d 


N(C 2 )/N(C 3 ) 






(xlO 12 cm- 2 ) 


(K) 


(K) 


(xlO 14 cm" 2 ) 


(K) 


(cm -3 ) 




HD 21483 




4.74 ± 0.28 
2.2 ± 0.5 f 


40 


560 


1.10 ± 0.30 
0.93 k 


13 ± 5 


190 ± 20 


23.2 ± 6.5 


HD 23180 


Per 


1.51± 0.13 
<1.55 f 


90 


470 


0.32 ± 0.12 
n 09 -t- n rn i 

\J.££ in u.uo 


60 ± 20 


200 ± 50 


21.2 ± 8.2 


HD 243Q8 


£ Per 


1 SQ _L n oc 

1.UJ _1_ U.UtJ 

1.0 s 


90 


235 


(I'll 4- n ns 

U.iJJ- _1_ U.UO 

n 1 q _|_ n no j 


80 ± 15 


ARC] + 1 c;n 

rlUU _1_ -LUU 


17.0 ± 5.8 


HD 24534 


X Per 


1.88 ± 0.19 
<X.5 * 


81 


408 


0.76 ± 0.20 
0.53 1 


44 ± 5 


225 ± 20 


40.4 ± 11.6 


HD 34078 


AE Aur 


6.5(3 ± 0.68 
2.2 ± 0.5 f 


171 


275 


1.00 ± 0.09 
0.58 k 


120 ± 10 


>500 


15.2 ± 2.3 


HD 36371 


v A nr 


0.43 ± 0.06 
<1.29 ' 


38 


554 


0.29 ± 0.11 






68.1 ± 27.9 


HD 53367 




2.21 ± 0.26 
<2.2 f 


40 


366 


0.61 ± 0.16 






27.6 ± 8.1 


HD 62542 




10.49 ± 0.53 


75 


137 


0.8 ± 0.2 11 


36 ± 15 


510 ± 50 


7.6 ± 1.9 


HD 148184 


X Oph 


2.50 ± 0.43 >> e 






0.43 ± 0.05 J 






17.2 ± 3.6 


HD 149757 


C Oph 


1.6 g 
1.50 ± 0.29 ! > e 






0.25 ± 0.07 
0.24 ± 0.03 J 






16.7 ± 5.5 


HD 152236 




1.8 ± 0.29 i ' e 


66 




0.16 ± 0.06 J 


36 ± 10 


300 ± 100 


13.3 ± 2.8 


HD 169454 




2.5 ± 0.66 
4.3 ± 0.5 f 


42 




1.60 ± 0.29 
0.70 ± 0.14 J 


50 ± 10 


>500 


64.0 ± 21.8 


HD 179406 


20 Aql 


2.0 g 
1.3 ± 0.7 f 






0.82 ± 0.15 
0.52 k 








HD 204827 




11.51 ± 0.87 
10.4 ± 0.5 f 


43 




4.40 ±0.29 


49 ± 5 


630 ± 200 


38.3 ± 4.1 


HD 210121 




3.79 h 
1.9 ± 0.5 f 






1.00 ± 0.13 
0.65 ± 0.15 m 









Ref eren ces. — (f 1 1Qka et alj 12003): fg) IMaier et al.H2 001D: fh ) IRoueff et aTl 120021): (i) IGalazu tdinov et all 120021) : HI Ivan Dishoeck fe Blacll 
1198611 : ( k) IFederman et al J 1199411 : f 11 iFederman fcTambertl ll98SD : fm) IGredelet all 1199211 : fn) IGredeT^HdTll993l) . 

Note. — Non-detections of C 3 in this work (with the upper limit JV(C 3 ) in 10 11 cm~ 2 in parentheses) arc: HD 11415 (1.8), HD 20041 (4.2), HD 
21398 (3.3), HD 22951 (2.5), HD 24912 (1.6), HD 351 49 (2.0), HD 370 22 (2.1), HD 37061 (3.3), HD 41117 (1.8), HD 42087 (2.1), HD 43384 (2.6), 
HD 206267 (3.5), HD 207198 (6.1), HD 210839 (6.1). IQka et al.| 12003H measure C 3 (with JV(C 3 ) in 10 12 cm~ 2 in parentheses) toward: HD 26571 
(2.1±0.6), HD 27778 (1.2±0.3), HD 29647 (4.6±1.3), HD 172028 (3.6±0.6), HD 203938 (<1.3), HD 206267 (2.7±0.4). 

a Temperature that best fits low J < 14 levels of C3. 

b Temperature th at best fits high J > 14 levels of C3. 

C AU values from Thorbur n et al.l 120031) unless otherwise noted 
and normalized to the oscillator strength / =1.0xl0 — 3 . 

d Temperatures and densities calculated from N\ values in Table 
□ 

c Reported equivalent widths have been converted to column 
densities using oscillator strengths in Table |5] and Equation 2. 
Weighted averages are used where multiple measurements are re- 
ported for the same J and the sum of the N j is reported here. Note 
that total N(C 3 ) reported here is roughly an order of ma gnitude 
larger than the values listed in Galazutdin ov et ail 120021) , which 
do not include the contributions of all observed rotational levels. 



levels of C3, with the kinetic temperature T{C2) deter- 
mined from the full excitation model of C2. While not 
unexpected, there is an excellent relationship between 
the kinetic temperature and the low J level population 
of C3. One would expect the temperature measured by 
fitting the higher ( J>14) rotational levels of C3 to be in- 
versely related to density, since as the density increases, 
collisional de-excitation is more efficient at reducing the 
high rotational level populations. This is indeed observed 
in the comparison of Thi g h{^3,) versus n co u and is shown 
in Figure 0] (right). 

It is interesting to note that the C3 columns measured 
toward HD 204827 and HD 169454, which are the two 
targets that have only low temperature components (Fig- 
ure |2l due to their high densities (n co /;>500 cm -3 ), are 
also the targets which deviate most from the trend in the 



relationship between Ti ow (Cs) and T(G2)- In these high 
density targets the kinetic temperature derived from C2 
and the temperature that best fits the excitation of C3 
are the same, whereas the low density (n co /;<500 cm~ 3 ) 
sources all show a somewhat increased excitation tem- 
perature of C3. This is because at high densities the 
excitation profiles of both C2 and C3 are predominantly 
determined by the kinetic temperature, whereas at low 
densities the low J populations, along with the high-lying 
levels, have a contribution from radiative pumping. 

5. SUMMARY 

We have observed rotationally resolved spectra of C3 
in 10 translucent sight lines. Using a new method of 
analysis to accurately retrieve the individual rotational 
level populations, comparisons have been made with the 
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Observations of Interstellar C2 Toward Various Stars 
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Line 


f.j',.j" 


HD 21483 a 


HD 23180 b 


HD 24398 c 


HD 24534 d 


HD 34078 a 


HD 62542° 


HD 152236 g 


HD 1694546 


HD 204827 h 


8757.69 


R(0) 


1.000 


9.3 ± 1.2 


1.3 ± 0.6 


0.57 


± 


0.17 


1.2 ± 0.2 




4 4-|_2 


0.7 ± 0.3 


0.7 ± 0.5 


5.6 ± 1 


12.8 ± 5.5 


8753.95 


R(2) 


0.400 


5.1 ± 1.6 


1.4 ± 0.6 


1.06 


± 


0.17 


1.9 ± 0.5 


4.2 ± 1 


9 ± 1 


1.6 ± 0.3 


1 ± 0.5 


6 ± 0.5 


18.3 ± 2.3 


8761.19 


Q(2) 


0.500 


11.5 ± 1.3 


1.5 ± 0.6 


1.19 


± 


0.17 


2.8 ± 0.2 


3.2 ± 0.7 


10 ± 1 


1.2 ± 0.3 


1.2 ± 0.5 


6.6 ± 0.8 


23.3 ± 1.8 


8766.03 


P(2) 


0.100 


5.6 ± 1.2 
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1 ± 0.5 
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1.4 ± 3 


8751.68 


R(4) 


0.333 






n so 
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OlOo. (0 




U.oUU 


*7 _1_ 1 Q 
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8750.85 


Rf6i 


0.308 




1.7 ± 0.6 


1.01 
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0.17 


3.2 ± 0.4 




2.8 ± 1 


1.7 ± 0.3 






7.2 ± 2.6 


8767.76 


Q(6) 


0.500 


8.1 ± 1.6 


1.8 ± 0.6 








2 ± 0.1 


4.6 ± 0.7 


1.9 ± 0.9 


1.4 ± 0.3 


1 ± 0.5 


3.9 ± 0.5 


14.3 ± 1.1 


8782.31 


P(6) 


0.192 












0.1 ± 0.3 












7.3 ± 1.7 


8751.49 


R(8) 


0.294 






0.37 
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0.17 










0.4 ± 0.4 






8773.22 


Q(8) 


0.500 




0.8 ± 0.6 












3.2 ± 0.9 


1.9 ± 0.3 




1.5 ± 0.5 




8792.65 


P(8) 
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0.286 
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8780.14 
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1 ± 0.8 


2.2 ± 0.8 
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Nj 


(xlO 14 ) 

















1.4 ± 0.2 


0.19 ± 0.09 


0.08 
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0.03 


0.17 ± 0.03 




0.6 ± 0.3 


0.1 ± 0.04 


0.1 ± 0.07 


0.8 ± 0.2 
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3.1 ± 0.3 


0.47 ± 0.14 


0.37 
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0.8 ± 0.06 
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2.9 ± 0.2 


0.45 ± 0.07 


0.36 ± 0.11 


2.1 ± 0.1 
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1 




2.1 ± 0.4 


0.5 ± 0.18 


0.42 


± 


0.04 


0.86 ± 0.1 


1.6 ± 0.2 


2.2 ± 0.3 


0.44 ± 0.09 


0.26 ± 0.15 


2.3 ± 0.2 


4.61 ± 0.7 




6 




2.4 ± 0.5 


0.61 ± 0.15 


0.48 


± 


0.08 


0.61 ± 0.04 


1.4 ± 0.2 


0.7 ± 0.2 


0.52 ± 0.08 


0.29 ± 0.15 


1.2 ± 0.2 


4.24 ± 0.3 




8 






0.24 ± 0.18 


0.19 


± 


0.09 


2.14 ± 0.13 




0.9 ± 0.3 


0.56 ± 0.09 


0.1 ± 0.1 


0.4 ± 0.2 






10 














0.35 ± 0.08 






0.21 ± 0.15 




0.3 ± 0.2 


0.64 ± 0.2 




n co ii 




190 ± 20 


200 ± 50 


460 


± 


150 


225 ± 20 


>500 


510 ± 50 


115 ± 20 


300 ± 100 


>500 


630 ± 200 




T i 




13 ± 5 


60 ± 20 


80 


± 


15 


44 ± 5 


120 ± 10 


36 ± 15 


43 ± 20 


36 ± 10 


50 ± 10 


49 ± 5 



References. — (al IFederman et al.l 11994ft : (Vi IHobbd 1198111 : fcl lChaffee et all 1198011 : (d) Private communication, Thorburn (2003); fei lGredel et all 1199311 : ff UHobbs fc CampbeTil 1198211 : (g) 
Ivan Dishoeck fc Blaclj I1989H : (hi This work. 
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Fig. 3, — Compariso n of the column densities of C3 and C2- N(Cs) determined in this work is plotted against N(C2) determinations of 
IThorburn et all 120031) along with best fit (dotted) line. Error bars are la. 
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Fig. 4. — Comparison of the excitation of C3 and C2. The excitation temperatures determined from the observed populations in low 
(J<14) levels of Cg are comp ared to the kinetic temperatures determined from the observed excitation of C2 and full excitation model of 
van Dishocck & Black (1986) (left), showing good agreement. The excitation temperature determined from the high (J>14) levels of C3 is 
compared to the densities derived from the C2 excitation model (right). As expected, Thigh (Qi) increases with decreasing density. 
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excitation of C2 , showing excellent agreement in the mea- 
surement of kinetic temperature. Trends with density are 
also consistent with the excitation expected for a non- 
polar linear molecule. 

With these measurements we have shown that the cor- 
relation between N(C2) and iV(C3), as pointed out by 
lOka et aD l)2003j) . is maintained for sight lines with a 
range of N(C 3 ) from 4.3xlO n cm~ 2 to 1.2x 10 13 cm~ 2 . 
The star HD 62542 does not fall on the correlation trend- 
line between N(C2) and JV(C3) and is shown to have 
either a very large column of C3 or depletion in C2, con- 
comitant with the known weakness of diffuse band fea- 
tures in its spectrum. Similar processes may be occurring 
to a lesser extent in AE Aur. 

Our measurements provide a new challenge for the fu ll 
C3 excitation model presented bv IRoueff et alJ {2002), 
which has been heretofore limited by the quality of ob- 
servational data. 
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